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ABSTRACT 

The predictabihty, or lack thereof, of the solar cycle is governed by numerous 
separate physical processes that act in unison in the interior of the Sun. Mag- 
netic flux transport and the finite time delay it introduces, specifically in the 
so-called Babcock-Leighton models of the solar cycle with spatially segregated 
source regions for the a and f2-effects, play a crucial rule in this predictabil- 
ity. Through dynamo simulations with such a model, we study the physical 
basis of solar cycle predictions by examining two contrasting regimes, one domi- 
nated by diffusive magnetic flux transport in the solar convection zone, the other 
dominated by advective flux transport by meridional circulation. Our analysis 
shows that diffusion plays an important role in flux transport, even when the 
solar cycle period is governed by the meridional flow speed. We further examine 
the persistence of memory of past cycles in the advection and diffusion domi- 
nated regimes through stochastically forced dynamo simulations. We find that 
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in the advection-dominated regime, this memory persists for up to three cycles, 
whereas in the diffusion-dominated regime, this memory persists for mainly one 
cycle. This indicates that solar cycle predictions based on these two different 
regimes would have to rely on fundamentally different inputs - which may be 
the cause of conflicting predictions. Our simulations also show that the observed 
solar cycle amplitude-period relationship arises more naturally in the diffusion 
dominated regime, thereby supporting those dynamo models in which diffusive 
flux transport plays a dominant role in the solar convection zone. 

Subject headings: Sun: activity — Sun: interior — Sun: magnetic fields 



Introduction 



Direct observ ations of su nspot numbers over 400 years, as well as proxy data for much 
longer timescales (iBeerl |2000| ). show that both the amplitude and the duration of the so- 
lar magnetic cycle vary from one cycle to the next. The importance of this phenomenon 
lies in the contribution of var ying levels of solar activi ty to long-term climate change, and 
to short-term space weather (INandy fc Martens! 120071 ). While there i s now a concensus 



that the Sun's mag netic field is generated by a hydromagnetic dynamo (lOssendrijveii 12003 



Charbonneaull2005l ) . the origin of fluctuations in the basic cycle is yet to be conclusivel y deter 
mined. Several different mechanisms have be e n proposed, including nonlinear effects ( Tobias 



1997; Beer et al. 1998 



Knobloch et al 



1998 



Kiiker et al. 



1999 



Wilmot- Smith et al. 



2005) 



2000 



Mininni fc Gomez 



Charbonneau et al 



2005 



stochastic forcing (lChoudhurilll992l : lHoyngill993l : lOssendriiver et al.lll996 



20021). and time-delay dyna mics (lYoshimura 



Char 



1978 



jonneau fc Dikpati 



Durney 



2000 



Wilmot-Smith et al.ll2006l ). A coupled, equally important, but 



ill-understood issue is how the memory of these fluctuations, whatever may be its origin, 
carries over from one cycle to another mediated via flux transport processes within the solar 
convection zone (SCZ). A unified understanding of all these disparate processes lays the 
physical foundation for the predictability (or lack-thereof) of future solar activity. These 
considerations motivate the current study. 

The main flux transport processes in the SCZ involve magnetic buoyancy (timescale on 
the order of months), meridional circulation, diffusion and downward flux-pumping (timescales 
relatively larger). Because magnetic buoyancy, i.e., the buoyant rise of magnetic flux tubes, 
acts on timescales much shorter than the solar cycle timescale, the fluctuations that it pro- 
duces are also short-lived in comparison. Our focus here is on longer-term fluctuations, on 
the order of the solar cycle period, that may lead to predictive capabilities. 
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Through an analysis of observational data, iHathaway et al.l (120031 ) have shown that the 
solar cycle amplitude and duration are correlated with the equatorward drift velocity of the 
sunspot belts during the cycle. They associate this drift velocity with the deep meridional 
counterflow that must e xist to balanc e the poleward flows that are observed at the surface 
( iHathawayl 1 19961 . l2005l : iMieschI l2005l ). The results show a significant negative correlation 
between the drift velocity and the cycle duration, so that the drift is faster in shorter cycles, 
consis tent with the interpretatio n of meridional circulation as the timeke eper of the solar 



cycle ( jNandyl l2004l : but see also ISchiissler fc SchmittI |2004| ) . In addition IHathaway et al. 



( 120031 ) identified positive correlations between the drift velocity of cycle n and the amplitudes 
of both cycles n and n + 2. While the two-cycle time lag was a new result, the positive 
correlation between circulation speed and amplitude of the same cyc l e is su pported by several 
earlier studies. In their surface flux transport model, IWang et al.l (120021 ) needed a varying 
meridional flow, faster in higher-amplitude cycles, to sustain regular revers als in the Sun 's 
polar field. They cited observational evidence from polar faculae counts (ISheeleyl Il99ll ). 
which peaked early for two of the stronger cycles, coinciding with poleward surges of magnetic 
flux. Furthermore, observations show a statistically-significant neg ative correlation between 
peak sunspot number and the dur ation of cycles 1 to 22 (Figure Ic of lCharbonneau fc Dikpati 
2000l : see also lSolanki et al.ll2002l ). Such a nega t ive co rrelat ion between cycle amplitude an d 
duration is also found in the models of iHoyngl (119931 ) and ICharbonneau &: Dikpatil (120001 ). 
Taken with the inverse relation between cycle duration and circulation speed, this is again 
suggestive of a positive correlation between circulation speed and cycle amplitude. 

Meridional circulation plays an important role in a certain class of theoretical solar 
cycle models often referred to as "fiux-transport" , "advection-dominated," or "circulation- 
dominated" dynamo models (see, e.g., the review by Nandy 2004). Such models have 



tures of the solar cycle ( 


Choudhuri et al. 


1995; 


Durnev 


1995 




Dikpati & Charbonneau 


1999; 


Kiiker et al. 


2001; 


Bonanno et al. 


2OO2I; 


Nandv & Choudhuri 


2001, 


2002 




Chatter jee et al. 



2OO4J ). In these models, a single-cell meridional circulation in each hemisphere (which is ob- 
served at the solar surface) is invoked to transport poloidal field, first poleward at near-surface 
layers and then down to the tachocline where toroidal field is generated. Subsequently, the 
return flow in the circulation advects this toroidal field belt equatorward through a region at 
the base of the SCZ which is characterized by low diffusivity. From this deep toroidal field 
belt, destabil ized flux tubes rise to the surface due to magnetic buoyancy, producing sunspots 
( lParkeiill955l ). We may point out here that the name "fiux-transport dynamo" is somewhat 
inappropriate to classify a circulation or advection-dominated dynamo (where the diffusion 
timescale is much larger than the circulation timescale throughout the dynamo domain). 
Our results indicate that diffusive fiux-transport in the SCZ could play a dominant role in 
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dynamos even when the cycle period is governed by meridional circulation speed, pointing 
out that flux-transport is a shared process. So, henceforth, by "flux-transport" dynamo, 
we imply a dynamo where the transport of magnetic field is shared by magnetic buoyancy, 
meridional circulation, and diffusion. 

Flux-transport dynamos offer the possibility of prediction because of their inherent 
memory. This arises specifically when the dynamo source regions for poloidal field production 
(the traditional a-effect) and toroidal field generation (the f2-effect) are spatially segregated. 
A brief discussion on important timescales (we identify three significant ones) in the dynamo 
process is merited here. The first is governed by the buoyant rise of toroidal flux tubes from 
the f2-effect layer to the a-effect layer to generate the poloidal field; since this is a fast 
process on the order of months, no significant memory is introduced here. The second 
involves the transport of poloidal field back into the f2-effect layer (either by circulation 
or diffusion). This could be a slow process where significant memory is introduced which 
is dominated by the fastest of the competing processes (advection versus diffusion). The 
third timescale relates to the slow equatorward transport of the toroidal field belt through 
the base of the SCZ, which sets the period of the sunspot cycle. In this class of dynamo 
models, with meridional circulation and low diffusivities in the tachocline (at the base of 
the SCZ), the third timescale is almost invariably determined by the circulation speed. It is 
the second timescale above, with competing effects of diffusive flux transport and advective 
flux transport, that becomes important in the context of t he persistence of memory . In th e 



advection-dominated, stochastically fluctuating model of ICharbonneau fc Dikpatil (l2000l ). 
this second timescale (governed by advection of poloidal field due to meridional circulation) 
was about 17 years, so that the polar field at the end of cycle n correlated strongest with 
the toroidal field of cycle n + 2 rather than that of cycle n + 1. The length of memory 
of any particular flux-transport dynamo model is unfortunately dependent on the internal 
meridional flow profile, and on other chosen properties of the convection zone which are not 
yet well-determined observationally. A particular problem is the strength of diffusivity in 
the convection zone, which strongly affects the mode of operation of the dynamo. 

Even if one assumes that these flux-transport dynamos capture enough of the realistic 
physics of the SCZ to make predictions of future solar activity, these predi ctions are critically 



depen dent on the relative role of diffusion and advection in the SCZ. iDikpati &: Gilman 



(120061 ). in their highly advection-dominated model, show that bands of latitudinal field from 
three previous cycles remain "lined up in the meridional circulation conveyor belt". They 
suggest that poloidal fields from cycles n — 3, n — 2, and n — 1 combine to produce the toroidal 
field of cycle n. Based on an assumed proxy for the solar poloidal fields (suns pot area), this 



leads them to predict that Cycle 24 will be a b out 5 0% stronger than Cycle 23 (IDikpati et al. 



20061 ). In stark contrast, IChoudhuri et al.l (120071 ). using a flux-transport dynamo model 
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with diffusion- dominated SCZ, and using as inputs the observed strength of the solar dipole 
moment (as a p roxy for the poloi d al fiel d), predict that Cycle 24 will be about 35% weaker 



than Cycle 23. IChoudhuri et al.l (120071 ) argue that the main contrib ution to the toroi dal 



field of cycle n, comes only from the polar field of cycle n — 1 (see also I Jiang et al.l 120071 for 
further details of this model). 

The confiicting predictions from these two solar dynamo models presumably result from 
the difference in the memory (i.e., survival) of past cycle fields in these models and could be to 



some extent infiuenced by the different inputs they use as proxies for t 



le so lar poloidal field. 



We also hypothesize that stronger diffusion in the IChoudhuri et al.l (120071 ) model destroys 
polar field faster, and that fiux transport by diffusion across the SCZ in this model short- 
circuits the meridional circulation conveyor belt, thereby shortening the memory of previous 
cycles. We perform a detailed analysis to test these ideas. To begin with, we consider a 
wider parameter space in the present paper, where we study the effect of varying meridional 
circulation speed and SCZ diffusivity on the amplitude and period of the solar cycle. In these 
simulations, we keep all other parameters the same, allowing a direct comparison between 
advection-dominated and diffusion-dominated SCZ regimes - which has previously been 
clouded by other differences between models. Then we introduce stochastic fluctuations in 
the model a-effect to self-consistently generate cycle-amplitude variations - as a completely 
theoretical construct towards studying cycle-to-cycle variations, in contrast to using diverse 
observed proxies for time-varying poloidal flelds. Subsequently, we perform a comparative 
analysis of the persistence of memory in this stochastically forced dynamo model in both the 
advective and diffusive flux-transport dominated regimes. Therefore, in spirit, this paper 
deals with the underlying physics of solar cycle predictability, and is not concerned with 
making a prediction itself. The layout of this paper is as follows. The main features of 
the model are summarised in Section [2l and the results of the parameter-space study are 
presented in Section [3J These results are interpreted in Section HI In Section [S] we analyze the 
persistence of memory in the advection versus diffusion dominated regimes. We conclude 
in Section [6] with a discussion on the relevance of this work in the context of developing 
predictive capabihties for the solar activity cycle. 



The Model 



We use th e solar dynamo code Surya , which has been studied extensively in different 



conte xts (e.g. iNandy fc Choudhuril |2002| . IChatterjee et al.l 12004 IChatterjee fc Choudhuri 



20061 ) ■ and is made available to the pubhc on request. The major ingredients of the code 
include an analytic flt to the helioseismically-determined differential rotation proflle, a single- 
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cell meridional circulation in the SCZ, different diffusivities for the toroidal and poloidal 
fields, a buoya ncy alKorithm to model radia l transport of magnetic flux, and a Babcock- 
Leighton (BL; lBabcocklll96ll . lLeightonlll969l ) type a-effect localized near the surface layer 
(signifying the generation of poloidal field due to the evolution of tilted bipolar sunspot pairs 
under surface flux transport). The code solves the kinematic mean-field dynamo equations 
for an axisymmetric magnetic field, which may be expressed in spherical coordinates (r, 9, 0) 
as 

B = 5(r,0)e^ + Bp, (1) 

where i?(r, 6) and Bp = V x [A{r, 6)e^] correspond to the toroi dal and poloid al components 
respectively. The mean-field MHD induction equation (see e.g. lMoffattlll978l ) then leads to 
the following standard equations for the a-Vt dynamo problem: 



^-^ + \{.-V){sA) = r^Jv'-^]A + aB, 



(2) 



(3) 



Here s = r sin 9, and we specify the meridional flow v, the internal angular velocity fi, the 
diffusivities rjp and rjt, and the coefficient a for the BL a-effect which describes the generation 
of poloidal field at the solar surface from the decay of bipolar sunspots. Note that although 
for modelling purposes the BL a-effect is mathematically similar to the traditional mean-field 
a-effect due to small-scale helical turbulence, the former is fundamentally different. The BL 
a-effect acts on much larger spatial (on the order of active regions or greater) and temporal 
(surface flux transport) scales, and is quenched at muc h stron ger field strengths (lO^G). The 
profiles of VL and a were described in IChatterjee et al.l ( l2004l ) and will not be repeated here. 
We will, however, describe the meridional circulation and diffusivity profiles in more detail. 



2.1. Meridional Circulation 

The meridional circulation is defined in terms of a streamfunction il){r,6), giving the 
velocity by 

pv = Vx[ij{r,e)e^], (4) 
where we assume the density stratification 



p = C {Rq/t - 0.95)^/^ . 



(5) 
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The streamfunction is given by 

/7]-(r - R 

ifjr sin 6 = ipQ {r — Rp) sin I — ^ 

\ Rq - Rp 

where /3i = 1.5 x 10'^m-\ p2 = 1-8 x IQ-^ m-\ e = 2.0000001, T = 3.47 x lO^m, and 
To = (-R0 — Rb)/4:. Here Rq = 6.96 x 10^ m is the radius of the Sun, R^ = 0.55/?© is 
the bottom of the simulation domain, and Rp = O.61i?0 is the penetration depth of the 
meridional circulation. We combine the arbitrary constants C and ipo in the parameter 
■^0 = — V^o/(0.95C) which gives, approximately, the flow speed near the surface at mid- 
latitudes. It is this parameter vq which we vary to change the circulation speed in this 
study. 

The circulation profile is illustrated in Figured] for vq = 25ms~^. The dots are plotted 
at yearly intervals for particles moving along the streamlines shown. 



2.2. Diffusion 



We use different diffusivities for the toroidal and poloidal fields, defined as follows: 



Vt[rj= Vrz + 



1 + erf 



BCZ 



+ 



1 + erf 



(7) 



Vp{r) = Vrz + 



V2 - Vrz 



1 + erf 



r - TBczj 



dt 



+ 



m - 112 



1 + erf 



rrcz^ 



dt 



Here dt = 0.025/?©, rscz = O.7-R0, r'^cz — O.72i?0, and ttcz = O.975i?0. In the radiative 
core we choose a low diffusivity, namely rijiz = 2.2 x lO^cm^s"^, representing molecular 
diffusivity only since there is no turbulent convection. We always choose rji < 772 so that the 
toroidal field diffusivity rjt in the convection zone is lower than the poloidal field diffusivity rjp. 
This is to model the suppression of turbulent diffusivity by strong magnetic fields, as toroidal 
field tend s to be strong and concentrated in localised flux tubes and therefore subject to less 
diffusion (IChoudhuril |2003| ) . whereas poloidal field is weaker and subject to more diffusion. 



At the surface both diffusivities increase to a high value (of the order of 10 cm s ), in 



line with surface flux transport models and observational estimates, 
illustrated in Figure [2j 



Typical profiles are 
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2.3. Numerical Domain and Boundary Conditions 

We solve Equations ([2]) and ([3]) in a meridional plane O.SSi?© < r < Rq, < 9 < 7r/2 
representing the Northern hemisphere. This is divided into a spherical gr id of 128 by 128 



cells, uniformly spaced in r and 6. We use the same boundary conditions as lChatterjee et al. 



( 120041 ). except that we consider only the Northern hemisphere and set i? = and dA/dO = 
at the equator (6* = 7r/2), thereby forcing the solution to have dipolar parity. 



2.4. Example Solutions 

Example solutions for two different runs are shown in time-latitude plots in Figures [3]^a) 
and (b), where the black lines denote contours of toroidal field strength at the base of the 
convection zone (r = 0.71Rq). This corresponds to the solar butterfly diagram, with the 
strongest field located at the active latitudes and migrating equatorward during each cycle. 
The background shading shows the strength of the radial field at the solar surface (r = Rq), 
which peaks at the pole several years after the toroidal field maxima (of the same sign) at low 
latitudes. The two solutions in Figure [3] characterize the diffusion-dominated SCZ (Figure 
I3t) and advection-dominated SCZ (Figure [3)d) regimes of the dynamo. In Figure MJo) the 
toroidal field shows a poleward branch at high latitude which is absent in Figure [3](a), and 
also a stronger radial polar field at the surface. The cause of these differences between the 
two regimes will become clear in Section HI 



3. Results 



We have carried out a parameter-space study to investigate how the cycle duration and 
amplitude in our model depend on the speed of meridional circulation and on the diffusivity 
in the convection zone. In each run of the code the parameters are held constant in time, 
but they are varied between different runs. Specifically, we vary the parameter vq, which 
gives the maximum circulation speed, and also ri2, which affects the diffusive decay and 
transport of the poloidal field in the convection zone, but not the toroidal field. In all runs 
we keep a surface diffusivity of ?7o = 2.0 x lO^^cm^s"^, and a toroidal field diffusivity of 
rji = 0.04 X 10^^ cm^ s~^ in the convection zone. These choices approximate the fact that 
turbulent diffusion is expected to be more efficient in the decay and dispersal of the weaker 
poloidal field, but less so for the stronger toroidal field; the latter suppresses the convective 
motions that give rise to turbulent diffusivity in the first place. The a-effect coefficient is 
not varied, but is set to ao = 30ms~^ for each run. This particular value was chosen to 
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ensure that periodic solutions could be obtained for a wide range of the parameters vq and 772- 
Each run was started from an arbitrary initial state, and then evolved until initial transients 
had disappeared, leaving a steady periodic dynamo solution. Such a periodic solution was 
found to exist only within a certain range of Vq for each value of 772- The cycle duration 
and amplitude were then measured from the periodic solutions, in the cases where such a 
solution was found. 

We define the cycle duration and amplitude by considering the time evolution of the 
toroidal magnetic flux $tor in a certain region around the base of the convection zone. 
Specifically, the toroidal field B{r, 6) is integrated at each time step over a region r = O.677i?0 
to O.726i?0, 6 = 45° to 80° (i.e., over the tachocline and latitudes 45° to 10°). This magnetic 
flux $tor should be proportional to the active region magnetic flux at the solar surface, 
under the assumption that more toroidal flux at the base of the convection zone leads to 
more buoyant eruptions. In a steady dynamo solution the flux $tor varies in strictly periodic 
manner, with its maximum amplitude giving the "cycle amplitude". We define the "cycle 
duration" to be the interval between successive peaks of |$tor| (half of the full dynamo 
period). This is therefore equivalent to the standard definition of the 11-year solar activity 
cycle, but of course the simulated periods are different. 

The resulting cycle duration and cycle amplitude are plotted as functions of the circu- 
lation speed Vo in Figures H] and [5] respectively. In these figures each curve corresponds to a 
different value of the diffusivity 772- The range of speeds covered by each curve indicates the 
range for which the code relaxed to a steady periodic dynamo solution, up to a maximum 
of Vq = 38ms~^. In Figure [6] the cycle amplitude is plotted as a function of 772, and in this 
case each curve corresponds to a different circulation speed vq. 



3.1. Dependence of Cycle Period on Meridional Circulation and Diffusion 



Figure H] shows a clear inverse dependence of the cycle duration on the meridional 
circulation speed Vq, with faster circulation leading to shorter cycles. A least-squares fit for 
the curve with 772 = 0.5 x 10^^ cm^ s~^ gives the dependence of the cycle period on meridional 
flow speed 

T = 217.716 ^;o"°•^^^ (9) 

found 



wh ere T is in years and f n is in rn etres per second. This agrees with the T ~ f g 



by iDikpati fc Charbonneaul (119991 ). and this inverse relation is a well-established result for 
Babcock-Leighton dynamo models. In these models the circulation, and specifically the 
equatorward counte rflow at the bottom of the convection zone, is the primary determinant 
of the cycle period (INandyl 120041 ). 
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The cycle duration is only weakly dependent on the diffusivity 772- A power-law fit for the 
curve with 772 = 2.0 x 10^^ cm^ gives T = 150.745 Vq^'"^^^ years. The lower power of fo here 
indicates that for higher- diffusivity solutions the cycle duration is slightly less dependent on 
circulation speed, presumably because flux transport by diffusive dispersal starts becoming 
important. Also, it is evident from Figure |4] that, at lower circulation speeds, there is a 
maximum diffusivity for which a periodic solution can exist. If there is too much diffusion 
at a low circulation speed, then the poloidal field will decay too much during its transport 
from high to low latitudes, thus generating insufficient toroidal field to sustain a periodic 
dynamo process. The essential difference between advective and diffusive flux transport is 
that the latter also reduces field strength during transport, due to diffusive decay. 

3.2. Dependence of Cycle Amplitude on Meridional Circulation and Diffusion 

Now we turn to the dependence of cycle amplitude on the speed of meridional circulation. 
This is shown in Figure El where each curve corresponds to a different diffusivity ri2 according 
to the legend. Rather than being monotonic, the cycle amplitude first increases with fo for 
low fo and then decreases with fo for large fo, with a turnover at some value of fo in between. 
The location of this turnover shifts to higher speeds as the diffusivity is increased. 

The dependence of cycle amplitude on diffusivity at any given circulation speed is not 
entirely clear on Figure [5l but is evident in Figure [6l where cycle amplitude is plotted against 
diffusivity 772. In this figure each curve corresponds to a different value of the circulation 
speed fQ. We see a similar behavior in that the cycle amplitude first increases with 172 for 
low ?72 and then decreases with rj2 for high ri2, with a turnover in between. If the circulation 
speed is increased, then the value of 772 corresponding to this turnover also increases. 

The behaviour of the cycle amplitude in our model, as illustrated in Figures [5] and 
[6l is more complex than expected a priori. Rather than a simple linear dependence on 
the circulation speed fo, there is a turnover in cycle amplitude, at a speed which changes 
depending on the diffusivity in the convection zone. In the next section we investigate the 
cause of this behaviour in the model. 



4. Advection versus Diffusion Dominated Solar Convection Zones 

The turnover of cycle amplitude as depicted in Figure occurs at a higher circulation 
speed fo as the diffusivity 772 is increased. The asterisks (joined by a thin line) in Figure 
[Tl^a) show the location of this turnover as a function of 772- We may think of this line in 
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the {'r]2,vo) plane as the dividing hne between two distinct regimes of the dynamo, which we 
call advection- dominated and diffusion- dominated. The advection-dominated regime corre- 
sponds to high circulation speed and low diffusivity, while the diffusion-dominated regime 
corresponds to low circulation speed and high diffusivity. A shift from one of these regimes 
to another affects flux-transport dynamics in a way that results in contrasting dependence of 
cycle amplitude on the governing parameters. Consider how the cycle amplitude varies with 
Vq for a fixed value of 772, corresponding to a curve on Figure O In the diffusion-dominated 
regime, a higher circulation speed means less time for diffusive decay of the poloidal field 
during its transport through the convection zone, leading to more generation of toroidal field 
and hence a higher cycle amplitude. In the advection-dominated regime, a higher circulation 
speed leads to a lower cycle amplitude because there is less time to amplify toroidal field in 
the tachocline (through which magnetic fields are swept through at a faster speed). It is the 
balance between these conflicting influences that leads to a turnover in cycle amplitude at 
some intermediate circulation speed. 

The bold line in Figure [Zl^a) shows the transition point between the two regimes that 
may be inferred from a simple balance between circulation and diffusion timescales tq and 
Tj). For a given circulation speed Vq, we define the circulation timescale tq as the time taken 
for meridional circulation to advect poloidal fields from r = O.95i?0, 6 = 45° to the location 
where the strongest toroidal field is formed at the tachocline [6 = 60°). The diffusion 
timescale is defined as td = L'^/ri2, where L = O.285i?0 is the radial distance across the 
convection zone from the same starting point. The two timescales are compared in Figure 
[7(b), where each horizontal line gives the circulation time tq for a different speed Vq, and 
the bold curve gives td as a function of 772. The crossing points of horizontal lines with this 
curve give the transition points between the advection dominated (rc < td) and diffusion 
dominated (rc > td) regimes from these simple theoretical considerations - which are in 
good agreement with the simulated transition points. 

4.1. Magnetic Field Evolution in Advection versus Diffusion Dominated 

Regimes 

We now compare the poloidal and toroidal field evolution in the two regimes. Figure [8] 
shows the poloidal field lines for two runs, at different times through the cycle, starting from 
one cycle minimum and finishing at the next cycle minimum, so that the fields reverse in sign. 
The left-hand column is taken from a run with vq = 20ms~^ and r/2 = 0.5 x lO^^cm^s^^, 
which is in the advection-dominated regime. The right-hand column is from a diffusion- 
dominated run with the same vo but with 7]2 = 2.0 x lO^^cm^s^^. Figure [9] shows the 
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evolution of the toroidal field for the same two runs. 

The key difference between the two regimes is the rate at which the poloidal field is able 
to diffuse through the convection zone, after it is generated at the surface by the Babcock- 
Leighton a-effect. This is seen clearly by comparing the poloidal field evolution between 8 
and 12 years in the two runs (Figures [Ht/h and d/i). In the diffusion-dominated run the new 
clockwise poloidal field diffuses directly down to the tachocline at all latitudes between these 
two times. However, in the advection-dominated run the new poloidal field does not reach 
the tachocline until the end of the cycle (16 years), and reaches the high latitudes before it 
reaches the tachocline; i.e., in this case, the field evolution follows the meridional circulation 
conveyor belt. There is still significant anticlockwise poloidal field remaining below the 
tachocline from the previous cycle, and even a lower band of clockwise field from the cycle 
before that. In the diffusion-dominated case there is only a weak band of anticlockwise field 
remaining from the previous cycle at solar minimum. 

This suggests that in the advection-dominated regime, poloidal fields from cycles n, n—1, 
and n—2 combine to produce the toroidal field for cycle n+1, while in the diffusion-dominated 
regime it is produced primarily from cycle n poloidal field, with a small contribution from 
cycle n — 1. 

The main difference in toroidal field evolution seen in Figure |9] is during the rising phase 
of the cycle, seen at 4 and 8 years (panels b/g and c/h). In the advection-dominated regime 
there are two separate regions of toroidal field generation, one region at high latitudes from 
poloidal field which has been advected poleward by the meridional circulation, and a second 
region at lower latitudes arising from direct diffusion of poloidal field across the convection 
zone. In the diffusion-dominated case there is no strong generation of toroidal field at high 
latitudes. In this case the strongest field generation occurs at mid to low latitudes, with direct 
diffusion presumably being the primary means of transporting poloidal field to the base of 
the convection zone (i.e., the meridional circulation conveyor belt is "short-circuited"). We 
point our however that, contrary to usual expectations, diffusive flux transport still plays 
a role in the advection-dominated case, and it is responsible for the complex dependence 
of cycle amplitude on diffusivity in the advection-dominated regime; we explore this issue 
below. 



4.2. The Role of Diffusive Flux Transport 

We have thus far identified the turnover in cycle amplitude to lie at the transition point 
between advection-dominated and diffusion-dominated regimes of the dynamo. Within the 
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umbrella of this model, this maximum in the amplitude is understood to be a balance 
between the time available for toroidal field amplification and the time available for poloidal 
field decay. Figure M shows how the cycle amplitude varies with the diffusivity 7]2 for a fixed 
circulation speed Vq. We see that there is a turnover in cycle amphtude for some value 
of T]2, with lower diffusivities corresponding to the advection-dominated regime, and higher 
diffusivities to the diffusion-dominated regime. In the diffusion-dominated regime, which is 
only reached for lower speeds Vq in Figure [6l the cycle amplitude decreases with increasing 
diffusivity. This is expected due to increased cancellation and decay of the poloidal field. 
However, in the advection-dominated regime, the cycle amplitude increases as the diffusivity 
ri2 is increased. This initially seems counter-intuitive, but we show here that it is caused by 
the influence of direct diffusive flux transport of poloidal field across the convection zone. 

This direct diffusion (from the solar surface to the base of the SCZ) was visible in 
the poloidal field evolution plots shown in the previous section; we now demonstrate its 
quantitative effect as 772 is varied, by comparing the poloidal field strength |Bp| at the base 
of the convection zone (r = O.TIS-Rq) with that at the solar surface (r = Rq). We take 
the ratio |Bp(base)| / |Bp(top)|, using the peak value of |Bp| at each location during the 
solar cycle. This ratio is plotted in Figure [TUl measured at latitudes 30° and 60°, and 
for two different circulation speeds. Thin lines correspond to Vq = 30 ms^^, where the 
dynamo is in the advection-dominated regime for the whole range of 772 shown. Thick lines 
correspond to vq = 20ms~^, for which the dynamo changes between the two regimes at 
about ri2 = 1.1 x lO^^cm^s"^ 

Consider first the behaviour at 30° latitude (solid lines in Figure [TOj) . Here the curves 
for both Vq have positive slope, which implies that a greater proportion of poloidal field 
from the surface reaches the bottom of the convection zone as ri2 is increased. Thus direct 
diffusive flux transport at lower latitudes always acts to increase the amount of poloidal field 
reaching the base of the convection zone. Nearer to the pole, at 60° latitude (dashed lines 
in Figure [To]), the behaviour is different. Here the ratio decreases as 172 is increased, both 
for the curve with vq = 30ms~^ and in the diffusion-dominated regime for vq = 20ms~^. 
In the advection-dominated regime for vq = 20ms~^ however, the ratio first increases with 
?72. This suggests a more complex relation between the surface and tachocline poloidal fields 
at high latitude, with competing influence from both diffusive and advective flux transport. 
This is expected because the downflow in the circulation is located at high latitudes. 

The analysis presented in this section supports the idea that direct diffusive transport of 
poloidal field across the convection zone, especially around mid-latitudes, is responsible for 
the trend of increasing cycle amplitude with increasing diffusivity, found in the advection- 
dominated regime. Although such diffusive transport acts to increase cycle amplitude in 
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both regimes, diffusion also causes the poloidal field that is being transported by meridional 
circulation to decay, cancelling with field from the previous cycle that is stored below the 
tachocline. Thus diffusion also has a negative effect on cycle amplitude. It is this negative ef- 
fect which dominates at higher diffusivities, forcing the dynamo into the diffusion-dominated 
regime where cycle amplitude decreases with increasing diffusivity. 



5. Persistence of Memory: Cycle-to-Cycle Correlations in Advection versus 
Diffusion Dominated Regimes in a Stochastically Forced Dynamo 

It is expected that the memory of a flux-transport dynamo is much longer in the 
advection-dominated regime than in the diffus ion-dominated regime, and solar cycle pre - 



dictions have been based on this expectation (jPikpati fc Gilmaru l2006l : I Jiang et al.l 120071 ). 
However, a detailed comparative analysis of persistence of memory in these different regimes 
under the umbrella of the same model had not been previously performed. Our analysis 
in the previous section has brought us closer to understanding the flux transport dynamics 
(in these two regimes) that is the physical basis for any memory mechanism. In this sec- 
tion, we consider how the persistence of this memory differs between the two regimes, by 
looking at the correlation between peak polar and toroidal fields of subsequent cycles. Since 
the simulations considered earlier relaxed to a regular periodic cycle, we cannot use these 
to study correlations between different cycles. Therefore, we now introduce self-consistent 
fluctuations in the cycle properties by means of a stochastically varying a-effect, and explore 
the resulting correlations between different cycles. 



5.1. Stochastic Fluctuations 



We introduce fluc tuations in the model by varying the coefficient ao of the a-effect (see 
Chatterjee et al.ll2004l for the full expression). We set 



ao = "base + «fluc 0-(t; Tc, 



(10) 



where abase = 30ms~^ is the mean value, ctfluc = 30ms~^ gives the maximum amplitude 
of the fluctuations (corresponding to the 200% level), and cr is a uniform random deviate 
selected from the interval [—1, 1], with a new value after each coherence time Tcor- Although 
for our purposes this is essentially a device for changing the cycle properties from one cycle 
to the next, there is a strong physi cal basis for stochastic variations in a which have been 



invoked in several previous studies (IChoudhuri 



Charbonneau &: Dikpatill2000l : lMininni fc G6meal2002l ) . Our model uses a Babcock-Leighton 



992 



Hoynglll993l : lOssendrijver et al.lll996 
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a-eSect where poloidal field is g enerated at the surface from the decay of tilted active regions 



flBabcocd 11961 



Leightonlll969l ). Thus stochastic variations in the a coefficient are natural, 



because it arises from the cumulative effect of a finite number of discrete fiux emergence 
events (active region eruptions with varying degrees of tilt). 

To compare the two regimes we consider two runs, both with ri2 = 1.0 x lO^^cm^s"^. 
The circulation speed Vq is kept constant throughout each run, and only the a effect is varied. 
Run 1 has Vq = 15ms~^, so is diffusion-dominated, while run 2 has Vq = 26ms~^ and is 
advection-dominated. The coherence time Tcor is set to 2.3 years in run 1 and 1.5 years in run 
2, so as to keep the ratio of the former to the cycle duration roughly the same in each case. 
We note that although the exact value of the coherence time is not important for our study 
(and is introduced just as a means to enable sufficient fiuctuations), the timescale - on the 
order of a year - is chosen to refiect that the BL a-effect is a result of surface fiux transport 
processes (diffusion, meridional circulation and differential rotation) which can take up to a 
year to generate a net radia l (component of the poloidal) field from multiple fiux emergence 
events (iMackav et aDl2004l ). 



5.2. Correlation Analysis 

In this section we compare the peak surface radial fiux for cycle n with the peak 
toroidal fiux $tor for cycles n, n + 1, n + 2, and n + 3. The toroidal fiux is defined as before 
by integrating B{r,9) over the region r = O.GTT-R© to O.726-R0, 9 = 45° to 80°. The radial 
fiux $r is found by integrating Br{RQ, 6) over the solar surface between 6* = 1° to 20°, (i.e., 
latitudes 70° to 89°). Note that the peak toroidal fiux precedes the peak surface radial fiux 
for the same cycle, which has the same sign. The poloidal field then produces the toroidal 
field for cycle n + 1 with the opposite sign. We measure the correlation of the surface radial 
fiux for cycle n with the toroidal fiux of different cycles, comparing the absolute value of 
each total signed flux. 

Both runs were computed for a total of 275 cycles with fluctuating ao, so as to produce 
meaningful statistics for each of the dynamo regimes. The results are illustrated in Figures 
[11] and [12] as scatter-plots of $tor for different cycles against $,.(n). The (non-parametric) 
Spearman's rank correlation coefficient is given above each plot, along with its significance 
level. The correlation coefficients are summarised in Table [T] where the Pearson's linear 
correlation coefficient is also given for comparison. Although the latter is less reliable, as 
it assumes a linear relation, it agrees well with rg in each case. 

The results show a clear difference between the two regimes. The advection-dominated 
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regime shows significant correlations at all 4 time delays, apparently suggesting that the 
memory of past poloidal field survives for at least 3 cycles; however, more on this later. The 
diffusion-dominated regime has a strong correlation only between ^r{n) and ^tor{n + 1), 
suggesting that the dominant memory relates to just a one cycle time-lag, although very 
weak correlations are also found with $tor(w) and $tor(^ + 3). 

In both regimes the strongest relation is the positive correlation between ^r{i^) and 
^tor(^ + This is to be expected as this is the more deterministic phase of the cycle 

- where toroidal fields (of cycle n + 1, say) are inducted from the older cycle n poloidal 
field via the relatively steady differential rotation. Note however that the two ffuxcs do not 
have to be directly coupled, in that the two fluxes may be positively correlated because 
they are both created from the mid-latitude poloidal field of cycle n (generated by the a- 
effect). The polar flux ^r{iT) arises through poleward meridional transport of the cycle n 
poloidal fleld, while the toroidal flux ^torin + 1) is generated from cycle n poloidal fleld 
that is diffusively transported across the convection zone. This is particularly the case in 
the diffusion-dominated regime. Nonetheless, even this indirect scenario suggests that the 
strongest correlation should in fact be between the cycle n poloidal fleld and cycle n + 1 
toroidal fleld in this class of a-fl dynamo models. 

The other phase of the cycle, in which the poloidal field is generated by the a-effect, is 
inherently more random due to the fiuctuating a-effect in these runs. Nevertheless, there is 
a strong positive correlation between $tor(^) and $,.(n) in the advcction-dominated regime, 
while this correlation is largely absent in the diffusion-dominated regime. This we attribute 
to the relatively stronger role of advective flux transport in the advection-dominatcd regime 

- which implies that a larger fraction of the original toroidal flux that has buoyantly erupted 
is transported to the polar regions by the circulation. In effect therefore, the advection- 
dominated regime allows correlations to propagate in both phases of the cycle, whereas the 
diffusion-dominated case allows correlations to propagate only in the poloidal-to-toroidal 
phase. The other correlation is broken in the diffusion-dominated regime because the advec- 
tion is short-circuited by direct diffusion, which transports flux downwards and equatorward 

- where it is cancelled by oppositely signed flux from the other hemisphere. This explains 
how the correlations can survive for multiple cycles in the advection-dominated regime, but 
not in the diffusion-dominated regime. 

6. Conclusion 

Signiflcant uncertainties remain in our understanding of the physics of the solar dynamo 
mechanism, implying that prediction of future solar activity based on physical models is a 
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challenging task. Here we have demonstrated how a flux-transport dynamo model behaves 
differently in advection and diffusion dominated regimes. Such differences, ar nongst others. 



have p reviously led to conflicting predictions of the amplitude of Cycle 24. iDikpati et al. 



( 120061 ) use an advection-domina ted model to predict a much stronger cycle than Cycle 23, 
whereas IChoudhuri et al.l (120071 ) use a diffusion-dominated model to predict a mu ch weaker 
cycle . The latter prediction is somewhat similar in spirit to the precursor methods (jSchatten 
20051 : ISvalgaard et al.l l2005l ) , which use the polar field at cycle minimum to predict the 
amplitude of the following cycle. Owing to the lack of observations of conditions inside the 
convection zo ne, opinions differ as to whether the real solar dynamo is weakly or strongly 



diffusive (e.g. iDikpati &: Gilmanll2006l : IJiang et al.ll2007l ) 



We find that for low circulation speeds f o (in the diffusion-domi nated regime), the cycl e 
amplitude is an increasing function of Vq, as in the observations of iHathaway et al.l (120031 ). 
However, the amplitude curve has a turnover point and is a decreasing function of vq at 
higher vo (in the advection-dominated regime), opposite to the observed correlation. When 
the diffusivity in the convection zone is increased, the location of this turnover moves to 
a higher vq. Our extensive analysis shows that this turnover corresponds to the transition 
between the diffusion and advection dominated regimes. In the diffusion-dominated regime, 
faster circulation means less time for decay of the poloidal field, leading to a higher cycle 
amplitude, whereas in the advection-dominated regime diffusive decay is less important and 
a faster circulation means less time to induct toroidal field, thus generati ng a lower cycle 



ampli tude. If the observed statistics of the past 12 cycles as reported by iHathaway et al. 



( 120031 ) reflect a true underlying trend, then our results imply that the solar dynamo is in 
fact working in a regime which is dominated by diffusive flux transport in the main body 
of the SCZ (although the cycle period is still governed by the slow meridional circulation 



counte rflow at the base of the SCZ). This conclusion supports the analysis of I Jiang et al. 
(l2007h . 



Through a correlation analysis in a stochastically forced version of our model, we have 
also explored the persistence of memory in the solar cycle for both the diffusion-dominated 
and advection-dominated regimes. It is this memory mechanism which is understood to lead 
to predictive capabilities in a-fl dynamo models with spatially segregated source regions 
for the a and Q effects. This understanding is based on the flnite time delay required for 
flux transport to communicate between these different source regions. We flnd that the 
polar fleld of cycle n correlates strongest with the amplitude (toroidal flux) of cycle n + 1 in 
both the regimes. In the diffusion-dominated regime this is the only signiflcant correlation, 
indicative of a one-cycle memory only. However, in the advection-dominated case, there are 
also signiflcant correlations with the amplitude of cycles n , n + 2, and n + 3. In contrast 
to the correlations that we infer, ICharbonneau &: Dikpatil (120001 ) found that the strongest 
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correlation in their advection dominated model was with a two-cycle time lag. Since such 
correlations lead to predictive capabilities, and obviously seem to be model and parameter- 
dependent as suggested by our results, such a correlation analysis should be the first step 
towards any predict ion, the latter being based on the former. In hindsight, however, both 
Dikpati et al.l (120061 ) - who use an adyectio n-dominated model and inputs from multiple 
previous cycles, and IChoudhuri et al.l (120071 ) - who use a diffusion-dominated model and 
input from only the past cycle to predict the next cycle, appear to be have made the correct 
choices within their modelling assumptions. 

Note that the memory mechanis m in our advectiori - domin ated case appears to have a 
different cause than that implied by iDikpati fc GilmanI (120061 ). who invoke the survival of 
multiple old-cycle polar fields feeding into a new cycle toroidal field. All of the surviving 
correlations in our advection-dominated regime (Figure [T^ are positive; they do not alternate 
in sign. This alternation in sign would be expected if bands of multiple older cycle poloidal 
field survive in the tachochne - odd and even cycle poloidal fields would obviously contribute 
oppositely because of their alternating signs. In that case we would expect the absolute value 
of $r(^) to correlate positively with $tor(^ + l) and ^torin + S) and s o on, but negatively with 
pterin + 2) and ^torin + 4) and so on, as evident in the results of ICharbonneau fc Dikpati 
(I2OOOI . Figure 9; after accounting for the fact that they use signed magnetic fields). Rather, 
in the advection-dominated regime of our model, the correlations appear to persist simply 
because fiuctuations in field strength are passed on in both the poloidal-to-toroidal and 
toroidal-to-poloidal phases of the cycle, as evide nced by the correlatio n betw een amplitude 
and polar fiux of cycle n. In a recent analysis, ICameron fc Schiisslerl (120071 ) find that the 
predictive skill of a surface fiux transpor t model - similar in spirit to the advection-dominated 
dynamo model of iDikpati et al.l (120061 ) - is contained in the input information of sunspot 
areas in the declining phase of the cycle. They argue that memory of multiple past cycles, in 
the form of surviving bands of poloidal field (its surface manifestations in their case), need 
not be th e only reason behind the predictive capability of the advection-dominated dynamo 
model of IDikpati et al. I J2006h. Our analysis of the advection-dominated regime supports 
this suggestion of ICameron fc Schiisslerl (120071 ). 



Coming back to the diffusion-dominated regime, our comparative analysis indicates 
that in this case, the memory of past cycles is governed by downward diffusion of poloidal 
field into the tachocline - which primarily results in a one-cycle memory. The fact that 
diffusion is an efficient means for transporting flux is often ignored, especially in this era of 
advection-dominated models; however, we flnd that diffusive flux transport is quite efficient. 
The identification of this one cycle memory in our stochastically forced model contradicts 
Dikpati fc GilmanI (120061 ) - who claim that prediction is not possible in this regime. As long 
as the source regions are spatially segregated, and one of the source effects is observable and 
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the other deterministic, flux transport a-Q, dynamos will inherently have predictive skills no 
matter what physical process (i.e., circulation, or diffusion, or downward flux-pumping) is 
invoked to couple the two regions. We may also point out tha t in the context of cycle-to- 



cycle correlations, downward flux pumping (ITobias et al.ll200ll ) would have the same effect 



as diffusion in that it also acts to short-circuit the meridional circulation conveyor belt. So 
although downward flux pumping differs from diffusive transport because in the latter case 
the fields may reduce in strength due to decay, the overall persistence of memory is expected 
to be similar if diffusive flux transport was replaced or complemented by downward flux 
pumping. 

In summary, our analysis has served both to explore the diffusion dominated and ad- 
vection dominated regimes within the framework of a BL type dynamo, and to demonstrate 
how the memory of the dynamo may be different in these two regimes. Based on our analysis 
we assert that diffusive flux transport in the SCZ plays an important role in flux transport 
dynamics, even if the dynamo cycle period is governed by the meridional flow speed. In 
fact, the observed solar cycle amplitude-period dependence may arise more naturally in the 
diffusion-dominated regime, as discussed earlier. Taken together therefore, we may conclude 
that diffusive flux transport is a significant physical process in the dynamo mechanism and 
this process leads primarily to a one-cycle memory which may form the physical basis for 
solar cycle predictions, if other physical mechanisms involved in the complete dynamo chain 
of events are well understood. Separate, detailed examinations of these other related physical 
mechanisms will be performed in the future. 
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Fig. 1. — Streamlines of the meridional circulation profile in our model (left), and latitudinal 
velocity profile across a radial cut at ^ = 45° (right). Dots on the streamlines show yearly 
positions for particles with vq — 25ms~^, starting from the squares at ^ = 45° and moving 
anti-clockwise. 
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Fig. 2. — Example diffusion profiles as a function of r for the toroidal (rjt) and poloidal (rip) 
fields. The dotted line shows the location of the tachocline. Here 770 — 2.0 x lO^^cm^s"^, 
r)i — 0.04 X lO^^cm^s"^ and r}2 = 1.0 x lO^^cm^s"^ 
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Fig. 3. — Two example solutions: (a) with Vq = 20ms~^ and 772 = 2.0 x lO^^cm^s^^ 
(characterizing a diffusive ffux-transport dominated SCZ); and (b) with the same Vq but 
r]2 = 0.5 X lO^^cm^s"^ (characterizing an advective flux-transport dominated SCZ). In each 
case black hues are contours of toroidal field B at the base of the convection zone (solid hues 
for positive values, dashed for negative). The grayscale in the background shows surface 
radial field strength Br{r — Rq), with white for positive and black for negative. The same 
contour levels are used in both plots. 



-26- 




Fig. 4. — Dependence of cycle duration on the meridional circulation speed vq. Each line 
style corresponds to a different value of 772 (the poloidal diffusivity in the convection zone) 
as given in the legend. 
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Fig. 5. — Dependence of cycle amplitude on the meridional circulation speed Vq. Each line 
style corresponds to a different value of 772 (the poloidal diffusivity in the convection zone) 
as given in the legend. 
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Fig. 6. — Dependence of cycle amplitude on the poloidal diffusivity 772 in the convection 
zone. Each line style corresponds to a different value of the meridional circulation speed vq, 
as given in the legend. 
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Fig. 7. — Transition between the advection-dominated and diffusion-dominated regimes. In 
(a) asterisks indicate the flow speeds Vq corresponding to turnover of cycle amplitude for fixed 
values of r)2 (inferred from the simulations shown in Figures 5 and 6). The bold line shows 
the transition point that may be inferred from simple theoretical comparison of circulation 
and diffusion timescales. Panel (b) shows the diffusion timescale td as a function of ri2 (bold 
line), and circulation timescales tc for selected speeds vq (horizontal lines), as defined in the 
text. 
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Fig. 8. — Comparison of poloidal fields in advection-dominated (left column/panels a to e) 
and diffusion-dominated (right column/panels f to j) regimes. Each row corresponds to a 
time through the solar cycle, running from one cycle minimum to the next. Solid lines show 
clockwise field hues and dashed hues show anti-clockwise field lines. Also indicated is the 
base of the SCZ at O.llRc,. 
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Fig. 9. — Comparison of toroidal field in advection-dominatcd (left column/panels a to 
e) and diffusion-dominated (right column/panels f to j) regimes. Each row corresponds 
to a different time through the solar cycle, running from one cycle minimum to the next. 
Grayscale contours show toroidal field strength, with black corresponding to the strongest 
negative field and white to the strongest positive toroidal field. Also indicated is the base of 
the SCZ at 0.71i?«. 
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Fig. 10. — Ratio of poloidal field |Bp| at base of convection zone (r = 0.715Rq) to that at 
the surface {r — Rq), measured as a function of diffusivity at latitudes 30° (solid hues) and 
60° (dashed lines). Thick lines correspond to runs with vq — 20ms~^ and thin hues to runs 
with Vq — 30ms~^. 
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Fig. 11. — Cycle-to-cycle correlations in the diffusion-dominated regime (run 1), between 
radial flux and (a) toroidal flux <l>tor(^), (b) ^torin + 1), (c) ^tor{n + 2), and (d) 

*^'tor(^ + 3). The Spearman's rank correlation coefficient is given along with its significance 
level for 275 cycles. All magnetic fluxes are in units of 10^^ Mx. 
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Fig. 12. — Cycle-to-cycle correlations in the advection-dominated regime (run 2), between 
radial flux and (a) toroidal flux <l>tor(^)) (b) ^torin + 1); (c) $tor(^ + 2), and (d) 

*^'tor(^ + 3). The Spearman's rank correlation coeflicient is given along with its significance 
level for 275 cycles. All magnetic fluxes are in units of 10^^ Mx. 
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Table 1. Cycle-to-cycle correlations 







(n) for run 


1 




(n) for run 


2 




(diffusion-dominated) 


( advection-dominated) 
















^tor(n) 


0.185 


99.8% 


0.287 


0.653 


100.0% 


0.778 


$tor(n + l) 


0.737 


100.0% 


0.706 


0.805 


100.0% 


0.851 


$tor(n + 2) 


-0.040 


49.1% 


0.028 


0.356 


100.0% 


0.546 


Pterin + 3) 


0.195 


99.9% 


0.202 


0.237 


100.0% 


0.417 


$tor(n + 4) 


0.036 


44.5% 


0.056 


0.183 


99.8% 


0.357 


^tor(n + 5) 


0.107 


92.3% 


0.013 


0.214 


100.0% 


0.316 



Note. — Correlation coefficients and significance levels for peak surface 
radial flux $r versus peak toroidal flux $tor for 275 cycles from stochas- 
tically forced dynamo simulations. 



